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Abstract 



'j^. An examination of the concept of using classical degrees of freedom to drive the evolution 

C^ ' of quantum computers is given. Specifically, when externally generated, coherent states of 

the electromagnetic field are used to drive transitions within the qubit system, a decoher- 
^ ' ence results due to the back reaction from the qubits onto the quantum field. We derive an 

expression for the decoherence rate for two cases, that of the single-qubit Walsh-Hadamard 
transform, and for an implementation of the controlled-NOT gate. We examine the impact 
of this decoherence mechanism on Grover's search algorithm, and on the proposals for use of 
error-correcting codes in quantum computation. 

PACS number(s): 03.67.Lx, 32.80.Qk 



1 Introduction 

In the original concept of quantum computation, the isolated, coherent evolution of a quantum 
system corresponded to a series of logical operations which could be used to compute a solution 
to a problem |l|]. In this concept, once a method of solving a given problem is decided upon, 
the logical steps in the method are translated into unitary transforms of a quantum system 0. 
These transforms demand a certain form for the Hamiltonian H that describes the evolution of 
the quantum computer, which in turn constrains the architecture of the quantum system. In other 
words, the program determines the system propagator, U{t) = exp{—iHt/h), which determines the 
form of H, whose parameters indicate the qubit-qubit interactions that must be present in order to 
carry out the program. 

Such quantum computers suffer from a significant drawback: they are not programmable. Con- 
sider the case of NMR quantum computing. A given program to solve a given problem determines 
H, which in turn determines the magnitudes of the J couplings between distinct spins, which 
in turn constrains the geometry of the molecule to be used as the quantum system. Once the 
molecule is created, perhaps utilizing organic synthetic chemistry, it is in general useful only for the 
method of solving the problem originally decided upon. In contrast with the flexibility of classical, 
transistor-based computers, these quantum computers have their programs "hard-wired" into their 
architecture. 

However, recent proposals for quantum computer architectures seek to overcome this limitation. 



They use external fields, generated by classical degrees of freedom, in order to drive the quantum 
system's evolution ||^. Since classical sources can easily be manipulated by the (presumably clas- 
sical) programmer, these methods offer a means by which the programmer can alter the evolution 
of the quantum system, and thus program the computer. Some of these proposals include the 
use of radio- frequency pulses acting upon nuclear spins in liquids [§|, ^ , laser pulses acting upon 
ions trapped in resonators [^, |^ , electrostatic fields generated from gated electrodes influencing the 
evolution of nuclear spins in semiconductors P], or electrons trapped in quantum dot structures 0. 
There is a question that naturally arises when one contemplates such proposals: how can the 
evolution of a quantum system remain coherent if it is interacting with classical degrees of freedom? 
Consider the loss of visibility of the interference pattern in a two-slit electron beam experiment 



which occurs when one attempts to measure which path the electron took using a light source [ITU 
The interaction of the field with the electron creates an entangled state of the electron and field, 
from which a measurement of the state of the light will reveal information about the position of 
the electron. To the extent to which the light carries information about which path the electron 



takes, interference is lost |T0|, [TI], 0. Such interference is central to the working of many quantum 



programs, e.g. Grover's search algorithm [0. If we use classically generated fields to drive (and 
thus interact) with the qubits, will that interaction effectively measure the state of the qubits, and 
by doing so, destroy the coherent evolution of the system? 

To answer this question, we examine a specific case, that of Grover's quantum search algorithm 



P^ . First, we review the implementation of this algorithm when quantum back reaction is not a 
concern. Then we assume that one of the transforms used in the algorithm, the Walsh-Hadamard 
transform, is driven by externally applied electromagnetic pulses. The quantum field is assumed to 
be generated by classical sources, so that it is described by a coherent field state. Thus, the quantum 
field exhibits behavior that is close to a classical field. Despite this, there is some quantum back- 
reaction on the field which leads to a decoherence of the qubit system. We determine the rate of 
decoherence due to this process. Finally, we discuss the broader implications of this mechanism of 



decoherence. While these questions have been previously raised [^, to the author's knowledge no 
quantitative assessment of their importance has yet been given. 

2 Grover's algorithm 

2.1 Discussion 

Before we consider how to drive a part of Grover's search algorithm, let us briefly review it here. 
Grover's quantum search algorithm is a method by which to retrieve elements in a subset of a larger 
set [|13|. It acts upon K qubits, or two-level systems, whose levels are arbitrarily label as |0) and 
|1). A complete, orthonormal basis for the system is the product basis, e.g. for ir=3, {|000), |001), 
|010), |011), |100), |101), |110), |111)}. It is usual to refer to each element of the basis set by the 
integer whose binary expansion represents the string of Os and Is, so that |3) = |011). The elements 



of the set to be searched are labeled by the integers from to 2^ — 1. 

The subset whose elements we are searching for is specified according to a condition. For 
example, we could search for any integer-valued roots of a given polynomials within a fixed range. 
There are many problems for which, given an x of the set of possible solutions, it can be checked in 
a polynomial number of steps whether x is a solution to the problem, but no known method exists 



to find all solutions in polynomial steps |15|. These problems can be solved by a brute- force search 
over all possible solutions, which is what Grover's algorithm does. Since Grover's algorithm has 
been shown to be optimal [|1^, 0, [1^ , its performance is one important indicator of how quantum 
computation might out-perform classical computation. 

In what follows, we assume there is only a single solution, y, to the problem. It is not dif- 
ficult to generalize this to multiple solutions [O. The initial state of the quantum computer is 
X]x=o ^ |x)/v2^. The goal is to transfer amplitude into the state \y) ||TB[ so that a measurement of 
the system yields the solution. This is achieved by a series of transforms of the form [WTZWO)^ . 
The unitary, Hermitian operator W = Hn^i''^) = Tin "\/2 {Sx{n) - Sz{n)) is the Welsh-Hadamard 
transform. It consists of a product of transforms, acting independently on each qubit in the system. 
(When referring to a single qubit, we employ the Pauli spin operator notation as commonly used in 



magnetic resonance [19]. In particular, Sa = |1)(1| and Sp = |0)(0| are projection operators, and 
5_|_ = |1)(0| and 5_ = |0)(1| are raising and lowering operators for a single qubit.) 

The operators TZ = —1 + 2|0)(0| and O = 1 - 2\y){y\ are diagonal in the product basis. The 



operator O is called the oracle. It is the only means by which the algorithm has knowledge of 
the solution. It would be the routine the flips the amplitude of the state \x), for example, if 
X were the root of a given polynomial. The combination WTZW can also be written as — 1 + 
2(Z]x=o^ \^))(J^'i=o'^{^\)/'^^ ^ which is called the invert- about- average step. 

The algorithm can be understood as a combination of two inversions, the first about \y) and 
the second about the state with an equal amplitude for all basis states. The two inversions result 



in a rotation, transferring amplitude into \y) [|l^, |T8l. For any normalized state of the computer, 
J2l=o^ ax\x), the application of the transformation WTZWO alters the state as follows: 
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which is a rotation of the probability amplitude between \y) and all other states, with siny9 ^ 2^~^'/^ 
for large K. When ay and Y.xj^^yO'x have the same sign, the amplitude for the state \y) increases 
with every iteration, and decreases otherwise. 

Note that IZ and C, unlike W, are not products of operators acting independently upon each 
qubit. They require qubit-qubit interactions to implement. To see this, write IZ = —1 + 2Y[n<Sf^{n) 
= — 1 + 2n„,[l/2 — ^^(rT,)]. When the product is expanded, terms such as Sz{n)Sz{m) appear, which 
indicate the need for qubit interactions. 
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2.2 Adding Classical Fields 

Consider how to drive the algorithm utihzing externally applied fields. The Walsh- Hadamard 
transforms W(n) can be driven one qubit at a time. Some proposals include methods by which 
qubit-qubit couplings, and thus TZ and C, could also be driven ||, 0. We assume here that only the 
VV(?2) are externally driven. Because of the nature of decoherence, it is reasonable to expect that 
if further transforms besides VV(n) are driven, the decoherence rate can only increase from what is 
derived below. 

Suppose that there exists a field / qubit coupling of the form KE{t)Sxi where E{t) is the 
electric field amplitude, and k is the field / qubit coupling. We assume the field uses only a single 
polarization for simplicity. If the qubits are magnetic dipoles, the form of the coupling is unchanged, 
with B{t) substituted for E{t). To give our programmer the greatest possible control, it is usually 
assumed that each qubit can be separately driven by the field. There are two ways to achieve this: 
spatial resolution (as is usually the case for lasers acting upon ions in traps) or frequency resolution 
(as is used in magnetic resonance). In either case, the Fourier components of the pulses acting upon 
separate qubits do not overlap. If spatial resolution is employed, then to prevent the pulses from 
overlapping, different directions of the beams are used. If frequency resolution is employed, then 
the pulses are centered at different frequencies in order to select a given transition. 



The entire qubit / field system is then described by a Hamiltonian of the form 

H/h = ^^ujnSzin) + Y^ Jn,mSz{n)S;,{m) + . . . + K^£;„(r„,t)5^(n) (2) 

n n,m n 

where the J terms, and higher-order spin couphngs if necessary, are present to implement the TZ and 
O transforms. The 'W{n) transform acting on qubit n at position r^ is driven by the field E^. Assume 
that the pulses have a square envelope with a center frequency Un, and that nEnifn^t) ^ Jn,m so 
that during a pulse we can ignore the J couplings (but see note PDI). Then implementing W for 



each qubit amounts to finding the propagator for a single field / qubit interaction acting for a short 
time, 

Hn/h = uJnSz{n) + KEn{rn) cos( u;„t )S,j;{n) (3) 

We want to adjust these parameters so that we arrive at the result VV(?2) = \/2{Sx{n) — Sz{n)). 
Suppose we write Win) = exp{—iiDnSz{n)t) Unit), which amounts to transforming to a rotating 
frame with the field. Remove the rapidly oscillating terms to arrive at 

d / K \ 

-K:Un = -i f (^n - iOn)Sz{n) + - -E„(f„) >Sa;(n) j Un (4) 

The propagator for a time- independent Hamiltonian of the form a ■ 5 is given by exp{—ia ■ St) = 
cos(at/2) - 2isin(at/2) a ■ S/a. In our case, a = { fi;i?„(r„)/2, 0, Un — COn }, so to arrive at the 
correct result we require the field to be detuned below the qubit, COn = ^n - ^/v^, with a field 
strength KEn{rn)/'^ = —^/V^ and a pulse duration of n/Q.. Then 

W{n) = tV2 exp (^-^7r^5,(n)) {S,{n) - 5,(n)) (5) 



This is the form we seek, except for an extra phase factor. Since yV{n) is apphed uniformly to each 
qubit in the system, if the spread of frequencies of the separate qubits is not too large, it will be a 
constant factor for the entire quantum computer, and can be ignored. Otherwise, suitable delays 
built into the system can be implemented to achieve the same purpose. 

3 Inclusion of Quantum Back Reaction 

3.1 Description of the Initial Field State. 

Before t = 0, when the computation starts, the programmer creates pulses of the field that propagate 
towards and drive the qubits at separate times tj. This is accomplished in a classical manner by 
turning on and off a classical current source, j{f, t), which interacts with the field through the vector 
potential, Hp{t) = — J dfA{'f,t)j{f,t). For simplicity, a single polarization for j is assumed. The 
current is classical in the sense that no quantum back reaction on the current source was included 
in the interaction Hamiltonian. 

Grover's algorithm requires a series of pulses, created at different times, call them ti, to drive 
each qubit n. Let the current distribution that drives qubit n be jn{r,t). It has significant Fourier 
components over a range of A; G K(n); as discussed previously, the K(n) are mutually disjoint sets. 
To create pulse i in the series, turn j„ on for a short time At at time tj. The field state is then 
transformed by exp{-iHp{jn{r,ti) )At/h) = TlfceKW ^d^k)- ^h^ ^d^) = exp(zag1" + z^ttj:) are 
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called displacement operators. They separately transform each mode of the field according to the 
complex- valued argument zj: oc At J j{f,t) exp{ik ■ r) dr. The new version of Grover's algorithm 
then has the form (W7^WO)^ E'=o ' k) nT""''' 0^'"' n^eKCn) D^ki^ki^^ti)) Kac), where |vac) 
indicates the vacuum state of the field. The W operators are now joint operators over the qubit 
and field. 

The displacement operators produce coherent field states that have many classical properties. 
A list of some of these is given in Ref. [^. We require the following properties: Dj:\z)Dj:{z) = 



1^, and D^ {z)Ojj:Dj:{z) = CLj: + z, and D^\z)Oj^ Dj:{z) = CLj:'^ + z*. The above properties lead to 
the following useful identity. Define the positive and negative frequency electric field operators 



SHf,t) = {S-{f,t)y = X: ^ ^^ % e^^-'^-"^ (6) 

The total electric field operator, S = S^ +S^ , obeys Maxwell's operator equations for a source-free 
region. Given a function / that can be represented by a Taylor expansion, then the following holds: 



/ (^+(f,t), £-{f,t)) exp{-zHp{j{r',U) )At/h) 
expi~zHpij{r',U) )At/h)f (^+(f,t) + E+(f,t), S-{r,t) + E-{f,t)) (7) 



Although / is still an operator, it contains the c-function E^{f,t), which is the electric field profile 
that would be expected classically at position r and time t, resulting from a current distribution 
j{f',ti) acting for a time At. 
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3.2 The Walsh-Hadamard Transform with a Quantum Field 

The goal of this section is to find a solution for a single qubit / field interaction, >V(n). As before, we 
assume that the pulses have square envelopes with center frequency a;„, and again let the detuning 
be given by Aa; = ujn — oJn- We also wish to distribute the displacement operators into Grover's 
algorithm in the following manner: 

npH(i«) n /^,K^,>,t.)) 7^^ m^)(i.-i) n ^fc(%>,t.-i)) o--- 

" [ fceK(n) J " [ fceK(n) J 

n >VH(i2) n D^,{z^,{nM))\nY[\w{n){t,) [] /^g(^,>,ti)) O e'|^)Kc) 

"■ [ fceK(n) J " [ fceK(n) J ^=0 

We have indicated the times at which the W{n){ti) interaction begins, e.g. the time at which the 
pulse arrives at the qubit. In commuting each D past 7?,, C, and Wirn) with m ^ n, the z^ pick 
up a phase of exp(— zcut). This corresponds to the free propagation of the pulse while the other 
calculations occur. D acting on qubit ra at a later time will not commute past W{n) at an earlier 
time. Classically, however, the pulse has not yet reached the qubit, so the expectation value of 
the £^ are zero. The non-zero commuter refiects the interaction of the qubit with the vacuum, 
including, e.g. spontaneous emission. Quantum systems are usually chosen explicitly to minimize 
this decoherence mechanism, so we assume that we can ignore it here. 

Start with the interaction Hamiltonian of Eq. (|^), except with the quantum field operators 
replacing the classical field. Transform to frames rotating with the field and qubit by assuming a 
solution of the form W{n) = exp{—iHQ{n)t/h)W where Ho{n)/h = J2k&K{n)'^0jj^'^(lj^ + uJnSz{n). 
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This results in an equation for W, 

jW{t) = -^1 (5+He-"*C(^"'«,i) +'5-He-^-"*^-(r-;,t)) W{t) (8) 

The modes over K(n) in the operator S^ slowly dephase during the interaction of the pulse with the 
qubit. This gives the effect of the pulse envelope on the qubit, but makes a solution somewhat diffi- 
cult. However, our interest is to find only the lowest order departures from classical behavior. Here 
is one way to do this: first, insert a unit factor into the propagator, exp{—iHQ{n)t / h) Y[K(n) ^k^k 
W{t) riKfra) -^fc I vac). The D to the right commutes to the front of the entire algorithm, and we 
solve for U = D'^WD instead. From Eq. (|^, this results in replacing S^ with £^ + E^ in Eq. (||), 
where E^ is the classical field at the qubit. 

Two things have changed. First, the appearance of the classical field profiles in Eq. (§) means 
that part of U alters the qubit evolution according to the classical field profile. Photons are still 
absorbed and emitted by the qubit, but the process occurs in such a way so that the coherent 
state of the field is not altered (which is a very "classical" behavior). Second, U still contains field 
operators, but now they act upon the vacuum state of the field with the D to the left. U will 
create one-photon states (and higher orders as well) which are orthogonal to the vacuum state. 
Because D is a unitary transform, the displaced one-photon states are orthogonal to the displaced 
vacuum states that describe the pulses. Thus, U can entangle the qubit state with a new field state 
orthogonal to the original field state. This is precisely the description of a decoherence mechanism. 

The idea that the classical field represents the lowest-order behavior suggests we should expand 
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out f/ as a series, J2j Uj, for which j corresponds to the power of the quantum field operators (or 
equivalently, the power of vh). This results in 



d 



-f/o = -^^ {E^S4n)e'-"' + i?-5_(n)e— "*) U^ 



dt 
for the classical behavior, and for j > 1, 

Ui incorporates the field operators once, so it represents the lowest order quantum effects. Each Uj 
has a solution of the form SaUajit) + SpUpjit) + S+U+jit) + S-U-j(t), 

|f/+,o = -^f e-*E+t/,,o |^-,o = -^^ e--*E„"^.,o (9) 

and for j > 1, 

|^+. = -^^ e^^"* (^n^^/^. + Cf/A.-i) jU^, = -^^e~^^"' {e-U^, + Cf/^.-i) (10) 

At this point we add in the classical field, whose envelope is a square pulse. Thus, let E^ = 

J^^T^'i>QT^^t^ wherc (j) is the phase, and 2E is the field intensity. Let C± = (d/dt)'^ ± iAuj{d/dt) + 
{k,E/2)'^ be linear differential operators. We can then rearrange Eq. (P) and ( p!0D to read 

/:„f/,,o = 0, C+Ufs,o = 0, £_f/+,o = 0, C+U-,o = (11) 

13 



and for j > 1, 



C U = - 
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The initial conditions are Ua^ = Upfl = 1, and U±fi = 0, and Uj = for all j > 1; and for the first 
derivatives, dUafl/dt = dUp^QJdt = 0, and dU±fi/dt = —i{K,/2)Ee^''^, and dUa,i/dt = dU^^i/dt = 0, 
and dU±^i/dt = —i{K,/2)S^{fn, 0), and dUj/dt = for all j > 2. The solutions for j = are 
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(13) 



where the tip angle is given hy 9 = JAlj'^ + (kE)"^ /2. This is a restatement of the usual expression 
for a Bloch vector influenced by a monochromatic field. 

The solutions for j = 1 are linear in the creation and annihilation operators. However, the field 
state they operate on is the vacuum state, so for the lowest order effect we will only require the 
solution for the creation operators. They are reasonably complex, so we simplify by setting Au 
= r2/v2, kE = —n/^/2, and the pulse length to n/Vl to reproduce the classical Walsh-Hadamard 
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transform. Then, to lowest order, and still in the doubly rotating frame, our transform is given by 
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The functions that give the Fourier components for the shape of the one-photon "back-reaction 
field" are given by 
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4V2(x" - 1) 
for which x = (a; — a;„)/r2, the offset of a field mode from the center frequency of the pulse, in units 

of the detuning. 

The system is assumed to be contained in a resonator of size L^ . The resonator is not perfect, 
but has a finite bandwidth associated with it due to coupling and conductive losses, which allows 
for the passage of the classical square pulse (which has the Fourier components (e*'^^ — l)/x ). 
Because all the g functions have significant Fourier components near Uj and decay as a;~^ like the 
square pulse, the resonator also allows the one-photon states created by the qubit to escape. The 
time-domain envelopes of these one-photon states (the Fourier transform of the g functions) are 
shown in Fig. [^. 

As in the classical field case, the extra phase factors of ±?7r/-\/8 in the transformation Eq. (|T^ 
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can be dropped, since they can be compensated for by proper timing. Similarly, we assume = 
at the start of the interaction. The (f) term does indicate a curious feature of the back reaction: g^ 
contains no term from the classical field, while gj^ requires two interactions with the classical field 
in order to perturb the coherent state. This is consistent with the property that the removal of a 
photon from a coherent state does not alter it, while the addition of a photon does. 
We can rewrite the transform Eq. (|14D in terms of unitary field operators. 
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(15) 



where the G operators create normalized one-photon field states, so that G'^G = 1. The normaliza- 



tion factors were found by numerical integration, for which the slowly varying ^/uj term was ignored: 
/^ = /^ \g^{x)\'^dx = 4.297, //3 = 4.297, /+ = 0.617, and /_ = 10.451. Although the G|vac) states 
are orthogonal to the initial coherent state, they are not mutually orthogonal. Later on, we will 



require their (non-normalized) overlap integrals Ia,i3 = -y//a//3(vac|G'J^G'^|vac) = J^^ g^{x)g/3{x)dx 
= 0.614 + i2.221, 4,+ = -0.617 + il.llO, /„,_ = 4.300 - i3.331, If3,+ = 0.617 + il.llO, J/3,_ = 
-4.300 -i3.331, and /+,_ = -1.850. 

The important result of this section is the transformation, Eq. (pISj). Each qubit / field inter- 
action has a probability amplitude, proportional to Q^^, to entangle the qubit with a field state 
orthogonal to the original field state. The greater Q, the lesser the decoherence via this mecha- 
nism. Actually, this may seem counterintuitive. Consider the interference pattern produced by a 
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coherent beam of electrons incident upon a double slit. Now allow a laser to interact with one of 



the two paths the electron could travel from the slit to the detector. If photons are scattered out of 
the coherent modes into vacuum states, then the visibility of the interference pattern is degraded, 



as expected |[T0|. If, however, only stimulated processes are important, then the visibility of the 
interference pattern should increase as the laser intensity is increased. The Poisson statistics of a 
coherent state can more efficiently hide the information about which path the electron takes as the 
number of photons in the beam increases. A similar situation has been noted with regard to welcher 
Weg experiments in atomic interferometry (see e.g. p2|). 



3.3 Grover's Algorithm with a Quantum Walsh-Hadamard Transform 



Taking the result from the previous section, Grover's algorithm is now 
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The displacement operators are not shown. They have been commuted all the way left, picking up 
the appropriate phase factors representing the propagation of the free field. The A and B are from 
transform Eq. (|15[), from which a factor of VV(n) is first removed: 
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The fact that operators acting independently upon separate qubits and field modes mutually com- 
mute was also used. 

When Eq. (|16|) is expanded out and only the lowest order terms are kept, to the original 
Grover's transform are added new terms. As previously discussed, the field operators in front of 
each term, one for each qubit and each pulse, create field states that are both mutually orthogonal, 
and orthogonal to the original field state. Thus, the probability that the final qubit state is \y) is the 
sum of the probabilities for each of separate terms to be in \y). In order to preserve normalization 
to (9(^2^^), the total probability to be in an orthogonal field state is subtracted from the probability 
of finding the field still in the original state at the end of the algorithm. The goal of this section is 
to calculate the success of this implementation of Grover's algorithm. 

After i successful steps of the algorithm, the computer state is given by 

cos(jV) v^ I \ ■ / ■ \ \ \\ \ \ 

I 2^ \x) + sm{jip) \y} |vac) 

V2^ - 1 x^y J 

where simp = 2-\/2^' — 1/2'^ [|l^. The general trend for the infiuence of the back reaction can be 
discerned from the specific example of 3 qubits and y = 2. A back reaction on the field can occur 
during either application of W. If it occurs immediately after O has been called (the operator B), 
or after the invert about average transform (the operator A). Suppose the reaction occurs after O, 
and it is for the least significant qubit. Then the computer state becomes 
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(5^5) + fi_|4)) + (5^16) + 5+|7)) + (5J7) + fi_|6)) - sinUv) ( Bp\2) + 5+|3) |vac) (1 



The amplitudes of pairs of states that differ at their least significant digit such as (0,1), (2,3), and 
so on, are mixed. If the error occurs for the second least significant qubit, then 



(5^10) + 5+12)) + (5,11) + B.\3)) + (5,13) + 5_|1)) + (5^|4) + 5+|6)) + 



cos(jV) 



(5,|5) + 5_|7)) + (5^|6) + 5+|4)) + (5,|7) + 5_|5)) ) - sin(jv9) (5^|2) + 5+|0)) Mvac) (19) 

The difference is in which pairs of states are mixed, and whether the qubit involved in the back 
reaction was initially in state or 1. Using this, we can sum over all the qubits of the computer, 
the total probability for a back reaction to occur during the first of the two applications of W, at 
step j in the algorithm. Suppose that the qubit frequencies are sufficiently close so that all pulse 
center frequencies are ui. Then 



iK-\\y\ 



cos(jV) 



K2^ -2 



2W^l ^°' (•^''^) ( l(^^ + 5_)|vac 



cos{jip) 



5_|vac) — sin(j(/9) 5^|vac) 



cos(jV) 



V2 



K 



5^|vac) — sm(j(p) 5_|vac) 



+ 






(B^ + B^Wsic) 



5, 1 vac) — sin(j(/9) 5+|vac) 



+ 



V2 



K 



5+ 1 vac) — sin(j(/9) 5,1 vac) 



The factor \\y\\ appears since the back reaction depends upon whether a qubit was initially in state 
or 1. Now suppose K is large enough so that 2^^ ^ 1. The following approximations are useful: 

X]J=i'^cos^(JV9) ~ So cos'^ xdx/if = rc/Aif, and I]J=i'^ sin^(jv9) ^ vr/4(^, and X^J^i*^ sin(jy9) cos(Jv?) 
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^ l/2(y9. Again, keeping only the largest terms in K, the probability that a back reaction occurs 
for any qubit, and at any step, during the first of the two W transforms is given by 



_2K/2 

mj eo 8 



y(|(i?, + i?+)|vac)|^ + |(5„ + i?+)|vac)|^) 



+ (K - |||/||)(|5^|vac)|' + |5+|vac)|') + ||y|| (|5„|vac)|' + |5„|vac)|') (20) 

If the back reaction occurs during the second W in an iteration, then j starts at 2 (the first invert- 
about-average step is always carried out), and the sign of the amplitude for \y) is positive. Taking 
the limit for large K, the cross-terms that depend upon the sign of ay drop out, and the final 
expression is the same as above except with the B operators replaced by those of A. The total 
probability of a back reaction anywhere in the algorithm is then the sum of these two. Starting 
with the operators of Eq. (|l3), the matrix elements are expressed in terms of the normalization 
and overlap of the Gjvac) states, e.g. \ {Bp + 5_) |vac) p + 1(5^ + B+) |vac) p = 2(/„ + Ip + I+ + I_) 
+ 43fJe(/^^_ -|- /a,+)- Plugging in, the final probability to end up entangled with an orthogonal field 
state is 

How much do the orthogonal field states contribute to the correct final answer? First, note that 
the matrix elements of A and B all have similar magnitudes. Thus, they equally mix the state \y) 
with the state connected to it by flipping the qubit that experiences the back reaction. Early in 
the algorithm when sin(j(y9) <^ 1, this does not increase the amplitude in ay significantly, since the 

20 



state that mixes with \y) has amphtude cos{j(p)/\'2^ <^ 1. At later times, however, the amphtude 
of \y) is near 1, so the back reaction decreases the probabihty to be in state \y) by roughly half. 

Thus, after j iterations of Grover's algorithm, a back reaction will cause the result to appear 
as if only ^ j/2 iterations had been performed. Further iterations will continue to transform these 
entangled states, but this does not necessarily improve the final result. Recall from Eq. (|l]) that 
amplitude is rotated into \y) only when the signs of ay and J2x^y k) ^^^ the same. For large K, the 
amplitude ay is ~ sm{j(f) [Jlj^Gp + -y7^G'_)|vac). The other qubit states are entangled with field 
states that are partly orthogonal to this state, but the amplitude that lies along the same direction in 



the Hilbert space of the field is, for large K, cos{j(p) {Ip + //3^+ + I^^p + /-,+) /yip + /- + 2'Ste{Ij3-) 
= (— 2.232 + iO. 448) cos{j(p). The real part has switched sign, and so further iterations will actually 
remove amplitude from \y). It should be clear that the back reaction essentially scrambles the 
memory of the quantum computer, from which further iterations will not, in general, recover the 
correct amplitude before the computation ends. 

4 Quantum Back Reaction for the CNOT Gate 

In this section, we briefly examine what happens when a classical field is used to drive a logic gate, 
such as TZ, that entangles together separate qubits. As before, let us focus on a single example, 
that of the controUed-NOT (CNOT) gate. CNOT is defined as a transform acting on two qubits 
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that flips the second qubit only if the first qubit is 1: |00) — ^ |00), |01) — > |01), |10) -^ |11), and 
1 11) —>■ 1 10). As for the Walsh- Hadamard transform, there are numerous ways to implement the 
CNOT gate [§]. Perhaps the most straight forward method is the use of a selective pulse. What we 
mean by this is the following. 

The two qubits to be transformed by CNOT have a Hamiltonian 

H/h = UiS, + uj2Tz + JS-X + Y, ^OLk^OLk + ^^{S^ + ^^) (^^ + ^') ■ (22) 

k 

—* — * 

By iS ■ T we mean the qubit operator SxTx+ <SyIy+ SzT^. Assume that uj2 > uji ':^ \uji — c<j2| 3> 
J. In this case, the transitions driven by the field occur at ^ u;i ± J/2 + J^/4(a;i — UJ2) and k, 
ijj2 ± J/2 — J^/4(co'i — ^^2), They appear as a "doublet of doublets" , that is, two sets of two lines each. 
Since H is nearly diagonal for small J, the highest frequency transition corresponds to |10) *^ |11). 
Thus, a CNOT gate can be implemented by selectively inverting this transition with a square pulse 
of duration T ^ ^-n j J . {T also depends on the matrix element for the transition.) The long length 
of the pulse keeps its bandwidth small enough so that it does not significantly overlap with the 
other transitions in the system. 

As before, first assume that the pulse is described by a coherent field state, and commute the 
displacement operator to the left. This replaces E^ with E^ + i^e^*"^* in Eq. (^2]). Second, split 
H = He + Hq, where Hq contains the electric (or magnetic) field operators terms. Suppose we 
find Uc(t), the propagator corresponding to He- Then the total propagator for the system can be 
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formally integrated as 



U{t) = Uc{t)~l f dt'Uc{t-t')HQUc{t') -^ f dt' f dt" Uc{t~t')HQUc{t' -t")HQUc{t") + - ■ ■ 
Ti Jo ft Jo Jo 

(23) 
Computing the value of the first integral in this series is the goal of this section. As before, this 
integral can be physically interpreted as follows: the qubit evolves under Uc, coherently evolving 
under the influence of a classical field until time t', when a scattering event occurs due to the 
quantum nature of the field. After the scattering, the qubit / field system evolves as before until 
final time t. Higher order corrections then correspond to multiple scatterings (or, multiple field 
operators). To find Uc, first assume a solution of the form exp(— i X^^Q-fc^flfc) exp{—iiD[Sz +Tz]t) 
Uci, and then drop the rapidly oscillating terms. Then diagonalize the resulting Hamiltonian, 



/ 



H, 



ci 



{uJi+uj2)/2-iu + J/4: kE/2 kE/2 

kE/2 {ijJi-uj2)/2~J/A J/2 

kE/2 J/2 {-uji+uj2)/2-J/A 





kE/2 
kE/2 



\ 



\ 



I 



kE/2 kEI2 -(cJi+cj2)/2+tu + J/4 

to find eigenvalues A„ and eigenvectors f„. Plugging in the Uc operator into Eq. (p3|) , and as before 



keeping only the terms linear in the creation operators, the lowest order quantum back reaction is 
(in the rotating frame) 






'ik-r 



2eaL^ 



O^k Y. VnviiS^ + I-)vlv^e' 



A„T 



n,m=l 



i{uj - OJ + Am - A„) 
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The real parts of the field amplitudes of the one-photon states entangled with different qubit states 
for a few of the matrix elements of U are shown in Fig. ^ The parameters used to generate the 
figure were ui = 20, uj2 = 21, and J = 0.4 for the qubits, and kE = 0.01, and T = 100.6 for the 
classical pulse profile. The classical result, Uc(T), transfers 97% of the amplitude between states 
|11) and |10). These parameters are not meant to represent any specific spectroscopic method. 

For on-resonant transitions, the matrix element peaks strongly at T oc E^^. Thus, as before, 
weaker pulses lead to greater decoherence. The matrix element for the transition |01) —>■ |00) is 
found to have the greatest magnitude. It may be surprising that transitions that are not being 
driven by the external field can also have large decoherence rates. Recall that this method of 
implementing the CNOT gate actually drives all the allowed transitions, but is only resonant with 
one of them. This decoherence process continues to grow because these other transitions that can 
lead to photon emission are weakly but continuously excited. 

In summary, several conclusions are of interest. First, the field states that become entangled with 
the two-qubit states are orthogonal to the initial field state, but they are not mutually orthogonal. 
In fact, their projections onto one another have essentially random phases. Computations that 
use the phases of the qubit state will be scrambled by this process. Secondly, because the CNOT 
gate involves entangling separate qubits, the scattering process also involves creating mutually 
entangled qubit / field states. It appears reasonable to suggest that this kind of field / multiple 
qubit entanglement will result when such logic gates are driven by external fields. 
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5 Discussion 

5.1 G rover's algorithm 

When unitary transforms are driven by externally generated coherent fields in the manner dis- 
cussed above, a decoherence mechanism exists that, with each applied pulse, tends to scramble the 
computer's memory. This decoherence mechanism is slightly different from the usual environmen- 
tally induced decoherence, in that it increases as the number of times the programmer attempts 
to manipulate the qubit system coherently. In the case of Grover's search algorithm where the 
Walsh-Hadamard transforms are externally driven, the degradation of the correct response scales 
as Q.2K2^/\hu/tQL''){K/^f 

A criticism of this analysis might be in the specific choice used to implement ^{n). Whatever 
method is chosen, the field / qubit propagators still hold, and some back reaction must exist (but 
see below). In general, the degradation should scale as the number of times a qubit transform is 
driven. For Grover's algorithm, if no error correction routines are implemented, then the amplitude 
VL of the field will have to increase exponentially with increasing K in order to keep the error below 
a fixed bound. Clearly, this is not a scalable way to implement Grover's algorithm. 

How important is this decoherence mechanism to the different proposed quantum computer 
schemes? Let us employ simple order of magnitude arguments, and ignore for the moment the 
implementation of error-correcting codes. The total probability for decoherence can also be written 
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as 0.2 [{huj) / [Li^ t^E"^ / 2)] K2^/'^. The pre-factor is the ratio of the energy of a single photon, to 
the total energy in the pulse (energy density times resonator volume), or the reciprocal of the total 
number of photons in the resonator used to create a pulse. 

First, examine the case for using lasers to drive single ions or atoms. A recent experimental 
demonstration of a logic gate using trapped ^Be^ ions as qubits |23| used 1 mW pulses of ~ 10^"^ 



s duration at 300 nm. This corresponds to 10^^ photons per pulse. The very small pre-factor will 
not pose a problem for computations involving a polynomial number of steps with increasing K, 
but for Grover's algorithm this mechanism limits the number of qubits to ~ 70. 

Next, suppose we examine the case for NMR quantum computing. First, let us address the 
question of utilizing the signal from a large number of independent quantum computers. Let us 
assume that our sample can be prepared in the ground state (all spins initially down), and let us 
ignore the interaction of the computers with one another through the field. We note that these 
assumptions are hardly justified for real world systems, but assuming a finite temperature for 



bulk quantum computing causes separate difficulties that have been addressed elsewhere ||2J. The 



application of Grover's algorithm results in a final state 

l^)= n ((l-7V2)|t/«) + 7EG'(a^)M^^))|vac) 
1 = 1 ^ 

where G creates orthogonal field states that contribute little amplitude to the correct solution. The 
total signal is the sum from all of the quantum computers in the sample, (^1^1 I^j ( |y^*^) (y^*-*! ) |^) = N 
(1 — 7^/2). The point is that at low temperatures, the macroscopic decoherence rate is multiplied by 
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the total number of independent quantum computers in the sample. Typically, NMR uses i/ ?^ 10^ 
Hz, or a photon energy of 10^^^ J. Pulses are 100 W for 10 /is, for a total energy of 10^'^ J. Thus, 
there are 10^^ photons per pulse. This limits Grover's algorithm for NMR to ~ 140 qubits if we can 
do NMR on a single spin system; but for a micro-mole (10^^) of quantum computers, we are limited 
to roughly 25 qubits. If electron spin resonance is substituted as the spectroscopic technique, then 
V ^ 10^° Hz, pulses are 1 kW for 10 ns, and thus require 10^^ photons, a negligible improvement 
over NMR. 

If, as is usually the case, ^2 <^ cjn, then the number of photons required to generate yV{n) is 
proportional to Q/(K'^uJn)- Thus, physical systems with small values oiujn or k are the most resistant 
to the above decoherence mechanism. Unfortunately, such systems have other limitations: if Un is 
small, then the temperature of the system is required to scale with increasing K in an unfortunate 
manner [21|, while k, is small, then the time required to drive the qubits increases, which in turn 



slows the entire computation down. It seems that quantum computers driven or programmed by 
externally applied fields will face significant design trade-offs. 

Will driving qubits by externally applied, static electric fields [§ offer any significant advantages? 
(Similar proposals might be envisioned for NMR by applying different static magnetic fields along 
different directions to individual spins.) In the Coulomb gauge, such longitudinal fields are not 



independent degrees of freedom 25 1. Instead, the total Hamiltonian of such a system must include 



the particles whose charges give rise to the static fields. Thus, for the case of qubits driven by 
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electrostatic fields from electrodes, decoherence results from the field operators (now with half 
integer spin statistics) for the Fermi levels of each electrode becoming entangled with the qubit 
operators. Finding the decoherence rate in this case is complicated by the extra structure for 
the electrons in the electrode, but it would be surprising if some decoherence mechanism was not 
present. 

5.2 Error Correcting Codes 

Finally, we briefiy examine some possible means to reduce the decoherence, and some difficulties 
they may encounter because of the unusual nature of this decoherence mechanism. First, consider 
the use of error- correcting codes. They are a general method by which to reduce decoherence 
p6| , |27| , provided that the decoherence is a statistically independent process amongst the qubits 



within each code- word. For logic gates applied bit- wise to code- words [^, these error correcting 
codes should be able to correct for this decoherence as for any other environmentally induced 
decoherence. The difficulty comes in when one considers how to actually implement the detection 
and correction scheme. Because of the requirement that errors within code-words be independent 
between qubits, we can not expect to have qubit / qubit couplings within each code-word. For the 
encoding schemes for which the authors are aware, this seems to require the "ancillary" bits |^ 
to be driven by external fields in such a way so that multiple qubits interact with the ancillary 
bit at once. Doing this with external fields, however, leads to entanglements between the qubits 
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within the code-word and the external fields. This in turn means that errors developed during the 
error correction scheme itself will not have the required statistical independence. However, in order 
to better understand how this decoherence mechanism might influence error correcting schemes, it 
would be helpful to examine a complete program, such as Grover's algorithm with a specific code 
implemented. 

Perhaps more hopefully, it appears that there are independent means by which storing infor- 
mation in multiple qubits could reduce this decoherence mechanism. Note that the appearance of 
terms such as \Bp + B-\ in Eq. (^) shows that destructive interference can lessen the probability of 
a back reaction. It is known to be possible to quench spontaneous emission in multi-level systems 
P8|| , which is a process quite similar to the decoherence mechanism considered here. Thus, it seems 
likely that a similar design could be used to remove the lowest order terms in the expansion of Eq. 
(P3|). Then, the decoherence rate increases as the inverse square of the number of photons in a 
pulse, which can be considered as a significant improvement. 

Conclusions 

Quantum computers that use external, classically generated electromagnetic fields to drive the 
evolution of the system undergo a decoherence induced by the quantum back reaction to those 
fields. The probability for the quantum system to be degraded increases as the total number of 
externally driven transforms, and inversely as the number of photons in the pulse. For algorithms in 
which no error-correcting codes are implemented, and for which the number of pulses is required to 
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increase exponentially as the problem size increases, the field intensity driving the system will also 
be required to increase exponentially in order to bound the degradation of the response. And for the 
implementation of error-correcting codes, the use of external fields for error detection and correction 
gives rise to a decoherence that does not have the property of being statistically independent between 
separate qubits in the code-word. 
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Quantum Back-Reacion Field Envelopes 



2.5 
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Figure 1: FIG. 1. A plot of the envelope of the field that describes the one-photon state emitted 
by a qubit whose state changes as indicated. For comparison, the envelope of the classical pulse is 
also shown. Time is in units of Q. 
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Figure 2: FIG 2. A plot of the real part of the frequency- dependent matrix elements for transitions 
(from top to bottom) |11) -^ |11), |11) -^ |10), |10) -^ |01), and |01) -^ |00). The lines at top 
represent the transitions in the two-qubit system driven by weak fields. 
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